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Abstract 



We face the well-known gyrokinetic problem, which arises in the description 
of the dynamics of a charged particle subject to fast gyration for the presence 
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\^ . of a strong electromagnetic field. The customary approach to gyrokinetic 
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of Darboux theorem, leads to potential complications and ambiguities due 

r^ . to the fact that canonical coordinates are field-related. Here we propose an 

-(— > ■ 

innovative formulation to construct gyrokinetic canonical variables based on 
the introduction of a new definition of canonical transformation. The new 

•i-H 

'^ . approach permits to shed light on this often controversial issue. 
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I. INTRODUCTION 

The "gyrokinetic problem" regards the description of the dynamics of a charged parti- 
cle in the presence of suitably "intense" electromagnetic (EM) fields realized by means of 
appropriate perturbative expansions for the its equations of motion. The expansions are 
usually performed with respect to the ratio e = ri/L << 1, where L and r^ are respectively 
a characteristic scale length of the EM fields and the velocity-dependent particle Larmor 
radius i"l = jf-, with Vlg = ^ the Larmor frequency and w the orthogonal component of a 
suitable particle velocity [see equation (47)]. The goal of gyrokinetic theory is to construct 
with prescribed accuracy in e the so called "gyrokinetic" or "guiding center variables", by 
means of a suitable "gyrokinetic" transformation, such that the equations of motion result 
independent of the gyrophase 0, being the angle of fast gyration, which characterizes the 
motion of a charged particle subject to the presence of a strong magnetic field. 

The first author who systematically investigated the gyrokinetic problem was probably 
Alfven^ who pointed out the existence of an adiabatic invariant, the magnetic moment /i, 
proportional to p^ the conjugate canonical momentum to 0, in the sense: 

|ln/i~0(e). (1) 

After subsequent work which dealt with direct construction methods of gyrokinetic 
variables^^^^, a significant step forward was made by Kruskal^^ who, first, established the 
consistency of the Alfven approach by proving, under suitable assumptions on the EM fields, 
that the magnetic moment can be constructed correct at any order n in e in such a way that, 
denoting M such a dynamical variable, it results an adiabatic invariant of order n, namely 
in the sense 

j\nMr^O{e") (2) 

M = /i + £/ii + ...+£>„. (3) 

A modern picture of the Hamiltonian formulation which makes easier the formulation of 
higher order perturbative theories, was given only later by Littlejohn^^ in terms of a non- 



canonical Lie-transform approach, adopting a suitable set of noncanonical variables. As a 
motivation to his noncanonical approach, Littlejohn^^^^^ pointed out what in his views was 
a critical point of purely canonical formulations such as previously developed Lie transform 
approaches^^'^^, namely the ambiguity in the separation of the unperturbed and perturbed 
contributions in the Hamiltonian due the presence of the vector potential A in the canonical 
momenta. He showed that this difficulty can be circumvented by making use of suitable 
non-canonical variables independent of A and which include the canonical pair {(f),p^). 

The possibility of constructing canonical gyrokinetic variables has relied, since, on only 
two methods due respectively to Littlejohn^^^^^ and Gardner^'^. 

The first approach, and probably the most popular in the literature^'^^^^ is the based 
on the use of Darboux theorem which allows, in principle, the construction of canonical 
variables for an arbitrary differential 1-form. The canonical 1-form expressed in terms of the 
canonical variables is then obtained by applying recursively the so-called "Darboux reduction 
algorithm" as pointed out by Littlejohn, which is obtained by a suitable combination of 
dynamical gauge and coordinate transformations. 

The second approach, due to Gardner^'^, was based on a mixed- variable generating func- 
tion formulation. This was used to construct the gyrokinetic canonical variables by means of 
a sequence of canonical transformations. In particular, an attempt to extend this approach 
to higher orders in e was later made by Weitzner^^. 

It should be stressed that in both approaches the canonical coordinates are field-related, 
namely they depend on the particular geometry of the magnetic field fiux lines. As a con- 
sequence the construction of canonical variables is achieved only under certain restrictions 
on the magnetic field (for the Lie transform approach see in particular 22,23^ p^j, example, 
White^^ assumed only small deviations from the axis-symmetric toroidal geometry, restrict- 
ing himself to weakly chaotic magnetic fields (assumption of regularity for the magnetic 
field), while Weitzner^^, based on his canonical approach, conjectured that the explicit con- 
struction of the gyrokinetic canonical variables might not be always possible. He claimed, 
indeed, that the magnetic moment might not result single- valued for locally chaotic magnetic 



fields, such as those occurring in non-symmetric MHD equihbria, for instance in Stellarators 
(assumption of "quasi-symmetric" magnetic field^^). An imphcation of his conjecture would 
be that for closed or ergodic orbits the magnetic moment would not be any more an adi- 
abatic invariant. This conclusion, if proven correct, on one side might indicate a potential 
breakdown of gyrokinetic theory itself and on the other might point out a difficulty intrinsic 
to the field-related choice of the gyrokinetic canonical coordinates. On the other hand, since 
no such assumption of regularity is required for the magnetic field in customary noncanoni- 
cal formulations, this conjecture raises also a potential contradiction between canonical and 
non-canonical formulations. An open issue is, therefore, the possibility of constructing more 
general, field-geometry independent canonical coordinates. 

The goal of this paper is to propose a new solution to this problem. Our approach is 
based on a new definition of canonical transformations, which we call generalized canonical 
transformations. These transformations make use of a suitable set of superabundant gy- 
rokinetic variables, which include in particular the canonical pair {(f),p^). Basic feature of 
the new variables is that they can be defined independently of the magnetic flux lines ge- 
ometry and do not require the use of the Darboux reduction algorithm. Basic consequences 
are that, on one hand, no such regularity or quasi-symmetry assumptions are required for 
the magnetic field, contrary to the conjecture of Weitzner, while, on the other hand, the 
magnetic moment in these variables results in all cases a single- valued function and therefore 
a physical observable. 

The paper is organized as follows: we preliminary reformulate Hamiltonian mechanics, 
extending it to the case of superabundant coordinates. Then, we face the gyrokinetic prob- 
lem, showing how a canonical transformation can be written in terms of superabundant 
gyrokinetic coordinates, which are shown to obey a suitable form of Hamilton modified 
variational principle. 



II. GENERALIZED CANONICAL TRANSFORMATIONS 

The purpose of this section is to extend the concept of canonical transformation. In fact, 
the customary definition (see for example^^'^^) may resuh, in some instances, too restrictive 
for actual applications. 

As is well known, the canonical transformation conventionally concerns Hamiltonian 
systems {x, iJ(x, t)}, natively represented in terms of canonical variables x=(q, p) and 
obeying the Hamilton equations 

x= [x,if (x, t)] (4) 

being H{^,t) the a suitably regular Hamiltonian function. In the phase space F, of dimen- 
sion 2g (where g is the degree of freedom), spanned by the vector x = (q, p), a canonical 
transformation is usually defined as a C*^"-* — diffeomorphism (with n > 2) 

7c : X ^ X (5) 

to a phase space V having the same dimension of the initial space F and satisfying the 
symplectic condition 

I = M ■ J ■ M'^ (6) 

is the symplectic (Poisson) matrix 



where M = |^ is the Jacobian matrix and J = 



of dimension 2g x 2g. 

It is well-known that canonical transformations can be defined also in the extended 
phase-space with dimension 2g + 2, spanned by the vector (q,gg+i = t,p,Pg+i) [see^^] and 
characterized by a set of superabundant variables (canonical extended variables) which sat- 
isfy Hamilton equations. Such extended transformations can be regarded as a particular case 
of what we shall define as generalized canonical transformation, i.e., a C*-"-* diffeomorphism 
(with n>2) 



7g : x^xg = XG(x,t,a), (7) 

(with a a real parameter) for which the transformed phase-space Tq can have larger dimen- 
sion of the initial phase space F and is characterized by superabundant variables, subject 
either to Hamilton or finite terms constraint equations. More precisely, first we assume 
dim(rG') = 2g' + k > dim(r) = 2g^ with g' > g and letting x^ of the form x^ = (z,u), 
with z = (zi, ...Z2g') and u = (-Ui, ...Uj?) . Then, we require that z and u obey respectively 
the extended Hamilton equations (for i = 1, ..., 2g') 

d d 

—Zi{:>^,t,a)= Yl 'Jij--g^K{^G,t,a), (8) 

j=l,2g' J 

(with J[, the canonical Poisson tensor of rank 2g' and assuming that the temporal variable 
t is left invariant by the transformation) and k constraint equations of the form 

/,(z,u,t) = 0, (9) 

with s = 1, fc and /s(z, u, t) real C*^^^ functions. 

The components of the transformed state x^ will be here denoted as superabundant 
canonical variables and the corresponding "equations of motion" (8), (9) generalized canon- 
ical equations. 

It is immediate to point out that generalized canonical equation can be set in variational 
form in terms of a suitable constrained form of modified Hamilton variational principle, just 
as the usual canonical equations. To provide a straightforward example, which will be also 
used in the sequel, let us consider the case of a charged point particle subject to an EM field 
and defined the following generalized canonical transformation 

x= (r,p) ^xg = (r,p,v). (10) 

One can prove that that the superabundant state xg is an extremal curve of the action 
functional 

5(xg)= rrft£(xG,r,t) (11) 

Jtx 



with the fundamental 1-form: 

C{xg, r ,t)dt = dr ■ p — dtH{r, p, t) — [dr — vrft] 
where the Hamiltonian is 



^ a' 

P —772V A 

C 



(12) 



^(r,P,t) = ^ 



1 r g ^2 



p-^A 

c 



+ g$ (13) 



(with {A, $} the EM potentials). It follows that the corresponding Euler-Lagrange (E-L) 
equations coincide with the canonical equations 

x=[x,//] (14) 

plus the non-holonomic constraint 

q 

p = -mv A. (15) 

c 

Therefore x^ defines a generalized canonical state for the Hamiltonian system. 

III. HAMILTONIAN GYROKINETIC THEORY IN SUPERABUNDANT 

VARIABLES 

In this section we intend to show that the previous Hamiltonian system, under a suit- 
able assumption of "strong" EM field, can be represented by means of a appropriate set 
of superabundant canonical variables which are gyrokinetic, i.e. the new Hamiltonian re- 
sults independent of the gyrophase when expressed in terms of them. The other basic 
feature is that the the corresponding canonical coordinates can always be chosen to be 
field- independent. To construct the new gyrokinetic variables we shall follow a two-step 
approach. 

The first step consists in constructing a particular set of hybrid^ superabundant gyroki- 



^In some cases it is convenient to adopt, for Hamiltonian or Lagrangian systems, a set of variables, 
which are not Hamiltonian or Lagrangian or Newtonian. In such circumstances we shall speak of 
hybrid variables. 
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netic variables, here denoted as pseudo-canonical x.' = {r' , pr' , (f)' , p^') . For definiteness let us 
require that the Hamiltonian function takes the form 



ff(r,p,i) = ^ 



p-^A 

ec . 



+ U (16) 

e 



where e is a real infinitesimal and assume that the EM potentials $, A are analytic functions 
of e and can be represented in the form 

oo 

$= 5:e^<l>,(r,t) (17) 

oo 

A= ^£^A,(r,t). (18) 

In validity of this assumption the construction of hybrid gyrokinetic variables is well known 
and has been achieved by several authors (see for example^'*^^^). In this case the Lagrangian 
expressed in terms of gyrokinetic variables (gyrokinetic Lagrangian) reads 

£(y', r', , t) =r ■^A*(r', u', f,', t)+ (19) 

^ + l) ^'B' - ^v'2 - U*{v\ u\ n\ t), (20) 

where the hybrid state y' is defined as follows 

y' = (r',«',/x',0'), (21) 

and, ignoring for sake of simplicity^ corrections of order 0{e^), one can directly prove that 
there results as a consequence 

V{r',u',^i',t)=u'h' + VE + eVj„ (22) 

where v^ = cE' x h'/B' is the electric drift velocity and v^ = ^ ^ 
I ^V'-B' + (u'h' + v^) ■ (-u'V'b' + Vv^) > is the diamagnetic drift velocity, both evaluated 



^In the sequel we shall omit higher order correction in e. 



at the guiding center position. Here the notations are standard. Thus b' = B{r',t)/B{r',t) 
while the primes denote quantities evaluated at the guiding center position r'. In particular, 
n' is the magnetic moment evaluated at the guiding center position. Moreover, {A*, $*} are 
the effective EM potentials, which at first order in e read 

A*(r', u', w', t) = A' + -—w' [I + 0{e)] , (23) 

$*(r',M',w;',t) = $'[l + 0(£)]. (24) 

To construct a set of superabundant variables, let us introduce the conjugate momenta 

Pr' = -?^ = -A* ^ mv'+^A', (25) 

in terms of which the gyrokinetic Lagrangian becomes 

£(x', r , , t) =r ■pr'+ p^' - K{x', t). (27) 

The superabundant state x' = {r' , p^' , (p' , p^') is here denoted as pseudo-canonical and 
K{^ ^t) is the corresponding Hamiltonian 



K{^,t) = -p^>n + ^ 



Pr' - -A' 

ec . 



+ -$*. (28) 



Manifestly, the transformation 

X = (r,pr) -^ x'= (r',pr', 0',P^') (29) 

is not canonical, even in the generalized sense previously indicated. 

The second step concerns the introduction of a further transformation to a new set of 
superabundant gyrokinetic variables. In particular let us consider the transformation 

7 : X = (r ,pr) ^ X'= (r',Pr' , 0', P<p' , v') (30) 

where v' = -u'b' + v^+ev^ is here considered as an independent variable. We intend to prove 
that the gyrokinetic state X' is canonical in the generalized sense defined above, namely 

9 



its components satisfy either Hamilton equations with respect to a suitable Hamiltonian 
function or finite-terms constraint conditions. To reach the proof, we initially notice that, 
by construction [see. Eq.(25)], the following finite-terms constraint is satisfied by the vector 



1 
m 



Pv 



ec . 



(31) 



Furthermore, using (r', u',p^i, (j)') the E-L equations for r' reads, omitting higher order terms 
in e 



Jt 



emc 



^Pr' + m ( u'Vh' + Vw'e + ^ V'A' ) ■ r + 



(32) 



+p^,V'n - m {u'W + V'v^) ■ [u'h' + vy - ^V'$* = 0, 



from one can prove that it follows 



(33) 



r = V . 



(34) 



As a consequence, there results 



-4pr' + - v'A' ■ v' + p..vn - ^v'$* = 0, 

at ec e 



which can be cast in Hamiltonian form with respect to K{x.',t) 



d (9 ^^, , , 



(35) 



(36) 



where the partial derivative d/dr' is defined keeping pr' as a constant. Analogously, from 
(31) and (34) one obtains the Hamilton equations 



d , d ^^, , , 1 

-r' = -— irx',t =- 
at opr' m 



Pr' --A' 

ec 



(37) 



Finally, in a similar way, it is immediate to prove that also p^, and 0' obey Hamilton equations 

d I d 



df 



-i^(x',t)=0. 



(38) 
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l.^'^±.K(At)--a'.. (39) 

Thus, given the constraint condition (31), which estabhshes a finite-terms equation for 
v', the remaining variables x.'= {r',pr',(f)',p^') satisfy the Hamilton equations, with re- 
spect to the Hamiltonian function i^(x',t). As a consequence, the superabundant state 
X'= (r',pr', 0',P(}!)', v') establishes a generalized canonical transformation [see for example 
(7)] in the transformed phase-space Tg- Furthermore X' results, by construction, a gyroki- 
netic state and therefore the gyrophase is ignorable for the generalized canonical equations 
(31, 36, 37, 38, 39). 

The crucial feature of these variables is that the canonical coordinates r' are manifestly 
independent of any particular magnetic field geometry, which implies that the vector r' can 
be represented in the form r' = r'(q',t), being q' = {q'l-, q2i Qz) arbitrary, field- independent, 
curvilinear coordinates, such as for instance orthogonal Cartesian coordinates. As a conse- 
quence, no restriction is placed on the magnetic field geometry for the definition of these 
canonical variables, contrary to previous formulations^'^'^^. A fundamental consequence is 
that, by construction, the magnetic moment /i can always be defined in such a way to 
be a single- valued function with respect with to any angle-like coordinates q' = (q'i, ^'2, Q's) 
and therefore results, in a suitable gauge, a physical observable. In particular, contrary to 
the conjecture of Weitzner^^, the definition of the magnetic moment results independent of 
the magnetic field topology and does not require the existence of a single family of nested 
magnetic surfaces (quasi-symmetric magnetic field) ^^. 

IV. CONSTRAINED HAMILTON MODIFIED VARIATIONAL PRINCIPLE 

A further key aspect of the present approach is that the previous generalized canonical 
equations (namely 31, 36, 37, 38, 39) are necessarily variational, as pointed out in section 
II. In fact, it is easy to show that they follow from a constrained form of Hamilton modified 
variational principle. One can prove that this is provided in terms of the following gyrokinetic 
Lagrangian: 
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£'(X', r,0,t)=r -pr' + 



-/C(x',t)- 



(40) 



r — V 



Pr 



C€ 



where /C(x',t) is the corresponding Hamiltonian 



''^^■'^-^ 



ec 



l^'sP^' + -/ 



(41) 



and by definition v' = u'W + v^ + ev^. 

The E-L equations corresponding to (40) coincide with the previous equations of motion 
(31, 36, 37, 38, 39) and hence provide a variational formulation for them. 

Finally, we emphasize that the constrained Lagrangian (40) defined above can also be 
obtained directly from the Lagrangian of a charged point particle in the presence of a strong 
EM field previously defined [see equation (12)]. This is achieved by means of the gyrokinetic 
transformation defined in terms of generalized canonical variables, i.e. of the form 



c' + eri + ... 



Pr =Pr' +^Pri + •• 

V = v' + evi + ... 



(42) 
(43) 
(44) 



with perturbations eri,epri,e\'i to be suitably defined. This transformation can be con- 
structed using standard techniques (see for example^^^^^) and can be identified at any order 
in e with a generalized canonical transformation. In particular, to the leading order in e one 
can prove 



Pr 
V 



r = r + ep , 

Pr = Pr' + ITT'W' + ^p'-V'A', 



(45) 



w 



where ep denotes the Larmor radius 



ep 



w'xb' 



(46) 
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The remaining notation is standard. Thus we require 



v' = M'b' + v^ + ev^, 



w'=w' {e[ cos 0' + e'g sin 0') , 

4''-arctJ<-m\. (47) 



^ mc,,/ 

P4>' ^/^ > 

where w' is a vector in the plane orthogonal to the magnetic flux line. Therefore, by means of 
the generalized canonical variables X'= (r' ,pr' , 4>' , p^' , V) gyrokinetic theory can be directly 
represented in canonical form without recurring to the use of hybrid gyrokinetic variables 
as considered by most of previous authors^^^^"^ 

V. FINAL REMARKS 

Based on the concept of generalized canonical transformation a new set of superabundant 
canonical variables has been defined for the set of gyrokinetic variables, which can be used 
to represent the Lagrangian of a charged point particle immersed in a strong EM field. 
These canonical variables, contrary to the customary ones previously considered in the 
literature^'^'^^'^^, do not depend on the magnetic field geometry and do not require subsidiary 
restrictions for the magnetic field. As a consequence, gyrokinetic theory can be developed 
in principle at any order in e, by means of gyrokinetic transformation expressed in terms 
of suitable generalized canonical variables. The present formalism for its straightforward 
simplicity and handy construction appears susceptible of interesting new applications both 
in plasma theory and mathematical physics. 
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